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, Confined optical phonons are discussed for a semiconductor nanowire of the Ge (Si) prototype 

■ on the basis of a theory developed some years ago. In the present work this theory is adapted to 

a non polar material and generalized to the case when the phonon dispersion law involves both 
I linear and quadratic terms in the wave vector. The treatment is considered along the lines of 

. a continuous medium model and leads to a system of coupled differential equations describing 

' oscillations of mixed nature. The nanowire is modelled in the form of an infinite circular cylinder 

and the solutions of the fundamental equations are found. We are thus led to a description of long 
wavelength optical phonons, which should show a closer agreement with experimental data and with 
calculations along atomistic models. The presented theory is applied to the calculation of optical 
phonons in a Ge nanowire. We have found the dispersion curves for various optical phonon modes. 
, We also normalize the modes and discuss the electron-phonon interaction within the deformation 

^ ■ potential approximation. 

PACS numbers: 63.22.+m, 63.20.Dj,63.20.Kr 
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I. INTRODUCTION 



The study of semiconductor nanowires is of current importance both for device applications and fundamental 
physics, a point that has been recognized since a certain time agOiA More recently, the fabrication of Si and Ge 
^ , nanowires by means of different techniques has been stimulated in view of their relevance in device nanotcchnol- 
Q . ogy i ^'^'^i^i^ The physical properties of these nanostructures are being studied both theoretically and experimentally, 
O [ while the fabrication of field-effect transistors on the basis of this type of systems has been also reported iSiSiifliiiiiSii^ 
In the present work we are interested in the phonon properties of non polar nanowires of the Ge and Si prototypes, 
I ■ and, particularly, in the optical phonons. It is well-known that phonon properties play an important role in the 
J> . considered systems because of their significance for the analysis of various physical processes, such as, charge and 
' thermal transport, optical transitions, etc. Moreover, phononic engineering has been invoked for device application 
ON : purposeSf^ making it relevant to have a better description of phonon modes. 

■ The analysis of phonon properties in quantum-wires already has a long history. By applying the long-wavelength 
approximation, based on continuum approaches, the polar optical phonons were discussed by different authors more 
than a decade agoi^*iSiiLi& (see also Ref. IT^for a more recent paper on the subject). Considering the nanowires made 
of non polar materials, we must remark the existence of a close related subject: phonons in carbon nanotubeSf^SiSi 
, an issue revealing outstanding possibilities for present day nanotechnology. 
d ' In the case of nanowires of the Ge or Si prototypes the calculation of optical phonons by means of long-wavelength 
, continuum approaches involves interest by itself, and also for applications. In the framework of the present paper 
I ' we shall support calculations of the optical phonon modes on a phenomenological theory, which was proposed more 
' O _ than a decade ago for the polar optical phonons in the same kind of nanostructuresi2SiSi24 Using this approach 
5 ■ polar optical phonons in quantum-wires were studied in Refs. It^IH . In the Ref. |^ the same theory, but adapted 
to non polar optical phonons, was applied to the case of carbon nanotubes. In all the above mentioned works it is 
explicitly assumed that the optical phonons display a dispersion law quadratic with the wave vector, which, in the 
bulk semiconductor case is of the form tj^ = ojq — /3^fc^ and, according to this long- wavelength approximation, should 
be considered near the F-point of the Brillouin Zone (BZ)iS^ In the present paper the theory is generalized in order 
^ to take account of linear terms in the phonon dispersion law, i.e., we assume a more general dispersion law having 
5^ [ the form uj^ — uj^ ~ P'k ~ (in the bulk case). Hence, we propose more general equations that in Ref. |^ now 
involving linear terms, and also we show how they can be solved. It is easy to realize that the additional linear term 
allows a better description of the optical phonons near the BZ center, in closer correspondence with both experimental 
evidence and atomistic calculations involving discretized microscopic models. In some of the studied cases this linear 
term becomes essential, and this situation is actually realized in nanostructures based on Si and Ge. 

Then, the paper is organized as follows. In Section II we present the fundamental equations, discuss their physical 
meaning and other details of the approach. In Section III we show how the equations can be solved in the case of a 
wire of circular cross-section and infinite length. Section IV is devoted to a physical discussion of the results obtained, 
while the case of a Ge nanowire is addressed. 
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II. FUNDAMENTAL EQUATIONS 

We assume that the relative displacement vector u satisfies the following equation: 

- ujI)u = pPjy"^ ■ u - pf3^\/ X V X M + ip{a ■ V)u + ipb{\/ ■ u) . (1) 

Equation |^ is analogous to that introduced in Refs. ll7ll8ll2H2^l23 . just considering that the electric potential is 
zero (see also Ref . l2l|) . However, new terms were added (two last terms at the r.h.s. of which modify the nature 
of the differential equation and lead to the linear terms in the phonon dispersion law. The parameter p is the reduced 
mass density of the two atoms (Ge) in the unit cell, ujq is the optical phonon frequency at the T point of the bulk 
semiconductor, and u{r, t) is considered to depend harmonically on the time (~ exp(— Parameters (3^ and (3t 
(with dimensions of velocity) are introduced to describe the quadratic (parabolic) phonon dispersion and used to fit 
the bulk semiconductor dispersion curve. They were already considered in previous works. Now we are introducing 
new parameters in the form of vectors a and 6, not considered in previous works and leading to linear terms in 
the phonon dispersion law. The latter vector parameters are also determined by a fitting of the bulk semiconductor 
dispersion law for the optical phonons. Equation (QJ represents a system of three coupled partial differential equations 
assumed to describe the confined polar optical phonons in each segment of the semiconductor heterostructure. The 
solutions must be matched at the interfaces by applying the boundary conditions for it, which we shall analyze below. 

It is important to notice that all the four terms at the r.h.s. of Eq. may be expressed as the divergence of a 
tensor cry , casted as 

cTy = p{P\ - 2/3|,)(V • u)6ij + p0^ {uj^t + Ui^j) + ip {b,Uj + a^Ui) . (2) 

This tensor could be compared to the stress tensor used in the theory of continuous elastic media, and actually its 
divergence provides a force density. However, in the framework of the present approach, it should be considered as a 
phenomenological quantity, which is introduced ad hoc in order to provide a description of dispersive optical phonons 
in close correspondence with what is seen both experimentally and in atomistic model calculations. In this spirit we 
do not attempt to explain the tensor cr.y in terms of the elastic properties of the system, and tensor aij should not 
be considered a stress tensor in the strict sense of this concept. The suffixes z, j = 1, 2, 3, denoting the cartesian 
components, should not be confused with the i = used as a factor. Obviously, in all the equations above we 

assume an isotropic model for the semiconductor. The short notation Ui^j denotes the partial derivative with respect 
to Xj, i. e., Ui^j — dui/dxj. 

Equations 1^ and (0) shall be now applied to the nanowire case, assuming it is an infinite cylinder of circular 
cross-section with radius rg. The use of cylindrical coordinates (r, 9, z) becomes a natural choice, and the equations 
must be written in terms of these coordinates, as well as the tensor aij. In the theory of elastic continuous media 
the standard boundary conditions applied to any interface between two different materials are: (1) continuity of the 
displacement vector it; (2) continuity of the force flux, given by a ■ TV, where N is the normal unit vector to this 
surface. Solutions showing divergences at any point of the region must also be disregarded and, in particular, we must 
avoid singular solutions in the asymptotic limits r 0, oo. 

It is convenient to introduce the auxiliary quantities F = V x it and A = V • it. Thus, after simple mathematical 
manipulations, Eq. is transformed into the following equations: 

[v2+gi]A = 0, (3) 
[V^+g2]r = * ax VA, (4) 

fjrp 

where 

oj^ + i{a + b)-V] , (5) 

w^ + ib-V] . (6) 

Notice that the latter quantities are actually differential operators in close correspondence with the starting equation. 
Another important issue involving Eq. I^J concerns its coupling to Eq. (j2Jl . All these features make the differences with 
respect to previous treatments fiii^*^ and lead to mathematical complications as a result of our aim of introducing 
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linear terms in the phonon dispersion law. When applying the equations to the case of an infinite cylinder of radius 
rg filled with an isotropic non polar semiconductor, we may make use of the translational and rotational symmetries 
of the system, assuming that the solutions must be proportional to exp[z(n0 + kz)], where n = 0, 1, 2, • • • and k 
plays the role of the wave vector. Moreover, vectors a and b shall be taken along the z-axis, provided that, with 
this choice, we are led to the appropriate linear dependence on k for the phonon frequencies. Under such conditions 
1l(t) become parameters linearly depending on the wave vector k and given by: = {loq — cj^ — (a + b)k)/ and 
g2 =(cj§-a;2-6fc)/^2. 

Another point of central importance concerns the boundary conditions that should be applied. In this work we shall 
consider the nanowire as "free standing" , so the appropriate boundary conditions involve zero force at the cylindrical 
surface of radius tq. After rewriting tensor in cylindrical coordinates, and realizing that the normal unit vector 
TV must be taken along the radial direction (TV = e^), we should require the components cr^r, crsr and azr to be equal 
to zero at r = rg. Starting from Eq. (0), tediously long but straightforward mathematical manipulations lead us to 
the following results: 
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III. SOLUTION OF THE DIFFERENTIAL EQUATIONS 

Assuming A(r, 9, z) — f{r) exp[i(n6' + kz)], Eq. © leads to 



d^f Idf f 

where ^ = Q^r and Qj^ ~ q\ — k^ . Eq. H1U|) is identified as Bessel equation of order n, and its linearly independent 
(LI) solutions are usually denoted by J„(^) and F„(^), Bessel functions of the first and second kind, respectivelyJS. In 
our case we shall apply just the function Jn(C)i convergent in the interior of the cylinder with < r < rg. 

Solutions of Eq. an inhomogeneous Helmholtz equation for the vector quantity F, are more difficult to find. 
In the first place we again shall consider that r(r, 0, z) ^ exp[^(n^^ + kz)]. We have found a particular solution of 
Eq. Q), which is given by 



Tp = Ae, X VA , with A = -^y^ — , (11) 

where is a unit vector along the z-axis. 

Then, we just have to find the solution of the homogeneous Helmholtz equation Th, and the general solution shall 
be given by F = Tp + Th- As discussed in Ref. I^B, page 1766, the two LI solutions of the homogeneous Helmholtz 
equation for a vector function F, using cylindrical coordinates, are casted as 



HI 



Vui X and Th2 = qTV2ez 



dz 



(12) 



where the functions Vi {i — 1, 2) satisfy the scalar homogeneous Helmholtz equations (V^ 



-4) 



0, with solutions 



analogous to those of A. After the general solutions for A and F are found, we can find the corresponding solution of 
Eq. Q for u, which shall be obtained from 
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In the foregoing calculations the involved mathematical steps are long and tedious, but here we shall report just 
final expressions. We write u — F{r) exp[z(n6' + kz)] where F = F^e^ + FqBq + F^e^ and e^., eg, and e^, are the unit 
vectors of the cylindrical coordinates. Then, we have 

Fr{r) = J^(Q^r)Ci + ^J,:(gTr)C2 + ^J„(QTr)C3, (14) 
Fe{r) = ^Jn{QLr)C^-^Jn{QTr)C2-J:,{QTr)C^, (15) 
Fz{r) = Jn{QLr)Ci + Jn{QTr)C2 , (16) 

WL 

where = — and 7 = a/{0^ — The "prime" in the Bessel function denotes the first derivative with 

respect to the function's argument and the constants C; (i = 1, 2, 3) should be determined after application of the 
boundary conditions. As it was already remarked in Section II, we consider a "free standing" wire, and the boundary 
conditions read as 

(Trr|ro=0, (T0r|ro=O, <J zr\ro — ^ i (17) 

indicating that the internal fiux of forces through the boundary surface is zero. Taking Eqs. 0, ©, and H17() 
into account we are lead to the following system of linear homogeneous equations 

3 

^G,,Q=0, (18) 

where the nine coefficients Gij are explicitly given in the Appendix. The dispersion relations, determining the 
frequencies of the various phonon modes as functions of the wave vector k, are obtained after solution of the secular 
equation 

det [Gy ] = , (19) 

while the components of the vector F are given as 

Fr{r)/C = SiJ:{QLr) + ^S2j;{QTr) + ^S3J.n{QTr), (20) 

Qt Qrr 

Fe{r)IC = ^S.UQlt) ^ ^S^-^UQtt) - SMQtt) , (21) 

Fz{r)/C = ~ SiMQlt) + 52 Jn(QTr) , (22) 
Ql 

where the functions Si {i ~ 1, 2, 3) are given in the Appendix. 

IV. NORMALIZATION OF THE PHONON STATES 

The determination of the constant C is made by normalization of the oscillation modes, which requires to transform 
the classical field u into a quantum-field operator. This is formally done in the form u —t ii = ua„. k, where a„^ ^ are 
second quantization bosonic operators. Operator an.k (a|i j.) annihilates (creates) a "n, k" optical phonon. We must 
also transform the classical kinetic energy for the oscillations to a quantum-mechanical operator by applying the rule 

Wkrn^]^Pj ii* ■ ud^r ^^pui^ I u*-ud^r^Hph, (23) 
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where 

Hph^jpuj'^J (^u^ -u + u-u^^ d^r^^puj'^ J u* ■ u d^r (^al ,^an ,k + ■ (24) 

In Eq. (|24(l Hph describes the free phonons hamihonian operator, while its hermitian character is ensured by 
construction. Requiring the latter expression to be given in standard form for phonons of the type "n , fc" (actually, 

we must have w ^^nk) Hph = huin ,k (ali fc^", + 5) ' '^^^ finally write the normalization constant C by 



Cr.,k = ^ —, T^j—, (25) 

where L is a normalization length (taken along the z axis) and Mn^k — J {P* i''')^r{r) + Fg{r)Fe{r) + 
F*{r)Fz{r))rdr/\Cn.k\'^- The explicit expressions for Mn.k involve rather complicated integrations of the Bessel 
functions and may be calculated numerically for each value of k and lu corresponding to the possible oscillation 
modes. 

For optical phonons in a non polar semiconductor the electron-phonon hamiltonian is given through the deformation 
potential approximation and, for an isotropic semiconductor, is proportional to V - it. Then, it is interesting to discuss 
this latter quantity, which provides an estimative of the interaction strength. 



V. DISCUSSION OF RESULTS 



With the aim of illustrating the foregoing theory, we consider a Ge nanowire. The Ge physical parameters are given 
in Table I. 



TABLE L Ge parameters 



material 


ao[10"*cm] 




LOO [cm ^] 


(3t [W^cm/s] 


Pl [lO^cm/s] 


a[10^**cm/s^] 


6[10'** cm/s^] 


Ge 


5.66 


5.32 


305 


3.48 


3.55 


-6.26 


3.59 



The parameters (3t, Pl, ol and h were determined by the authors by a fitting of the bulk semiconductor phonon 
dispersion law reported in Ref^Jj considering dispersion along the (001) direction and embracing approximately 
50 — 40 % of the BZ near the F point. Parameter luq was taken from the same reference, while the other parameters 
can be found in Ref. l28l 

In Fig. 1 we depict the dispersion relation for several axis-symmetric modes (n — 0) taking rg = 10^'' cm. The 
frequency is in units of cm~^) and the wave vector k is in units of 27r/ao, where oq is the lattice parameter for the 
bulk semiconductor). The curves were obtained from Eq. H19|l taking n = 0. Let us emphasize that these modes 
do not present torsional oscillations, and they result from a combination of radial breathing vibrations together 
with vibrations along the cylinder axis (axial oscillations). The latter oscillations are not purely longitudinal as it 
is sometimes argued, but a mixture of LO and TO oscillations. However, at A: = the axis-symmetric oscillations 
separate into pure LO and TO vibrations, a fact we have used for the labelling of the different modes. 

Notice that for the curves corresponding to the modes TO2 and LO3 there is an interval of k values (approximately 
between 0.5 and 0.63) where the curves are interrupted, and no frequencies are reported. This region obviously 
corresponds to an anticrossing of the mentioned modes. Of course, in the real case, there must be oscillation frequencies 
and the interrupted curves may be interpreted as a limitation of our approach, which is not able to describe such 
oscillation modes. However, it is worth to mention that the observed failure has nothing to due with the anticrossing 
effect by itself. The theory we are applying is essentially valid not too far from the fc = region, and we are actually 
extrapolating it beyond this region. This latter issue should be seen as a singular failure of the applied treatment, 
which, otherwise is giving a rather good description of the dispersion law for these modes. We have set rg = 10^^ 
cm, a relatively small radius which, however, is typical for this kind of nanowires. For such thin nanowire the number 
of physically possible optical phonon modes is rather limited and, in fact. Fig. 1 displays a few more modes than 
actually present. In the framework of the current paper one important contribution is to stress the importance of 
taking the role of the linear term into account we have added to our equations. It is easy to convince ourselves of this 
latter issue: making the linear terms equal to zero the general structure of the dispersion law depicted in Fig. 1 is 
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wave vector {In/a^ 

FIG. 1: Dispersion curves for the axis-symmetric modes (n=0). The frequency is in units of cm~^ and the wave vector k in 

units of 27r/ao) and we take ro = 10^^ cm. 



300 




wave vector {In/a^) 

FIG. 2: Dispersion curves for the flexural modes for n = 1 taking ro = 10~^ cm. The frequency is in units of cm~^) and the 
wave vector k in units of 27r/oo). 

subjected to significant changes. The other interesting point is to compare with experimental results or calculations 
based upon microscopic (atomistic) approaches. Unfortunately, in the revised literature we have not found reports 
on optical phonons for this particular type of Ge nanowires. 

In Fig. 2 wc show the dispersion curves for the so-called flexural modes taking again ro = 10^^ cm and for ri = 1. 
These modes involve all kinds of possible vibrations: radial, axial and torsional. Torsional modes can be considered 
for any n > and we are analyzing the case with n = 1. In the chosen frequency interval (140 < co < 300 cm~^), 
involving the higher frequency values, just nine modes arc present. We should remark that, for the flexural modes, 
the oscillations have always a mixed character and cannot be separated into pure LO and TO oscillations. In contrast 
with the axis-symmetric modes the mixing is present even at fc = 0, and then we cannot classify the modes in the 
same way as we did in Fig. 1. Other feature clearly seen in Fig. 2 is the presence of several anticrossing regions: curves 
belonging to the same symmetry cannot cross each other and are "repelled" whenever they get too close. Finally, 
there are several curves (those with the higher frequencies) that are interrupted for the higher values of k (especially. 
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FIG. 3: Absolute value of the divergence of vector u (in units of <I>o) as a function of r/ro for the axis-symmetric modes (n = 0). 
We fixed k = 0.2(27r/ao) and considered the frequencies 174.8, 201.9 and 238.1 cm~^ corresponding to three possible modes 
displayed in Fig. 1. 




FIG. 4: Absolute value of the divergence of vector u (in units of $o) as a function of r/ro for the flexural modes (n = 1). 
We fixed k = 0.2(27r/ao) and considered the frequencies 173.18, 205.8 and 229.4 cm^^ corresponding to three possible modes 
displayed in Fig. 2. 



in the region k > 0.6, in units of 27r/ao). The explanation for this effect is the same as in Fig. 1: a failure of the used 
approximation, which is actually valid for values of k not too far from fc = 0. However, we must again remark that the 
theory actually gives a very useful description of the considered modes for a relatively wide interval of k (60 — 70 % 
of the BZ). 

Besides the dispersion curves, we find interesting to analyze the quantity given by V • it = div {u), as a function of r 
in the interval r < rg. The point is that this quantity is proportional to the electron-phonon interaction according to 
the deformation potential approximation, and then can give us an estimative of the strength of this interaction. We 
have taken the divergence of vector u, given in Eq. (|22|l in cylindrical coordinates, together with the normalization 
constant, Eq. H25|) . which was evaluated numerically. 

In Fig. 3 we plot the absolute value of the mentioned quantity {\div {u)\) in units of $o = [Ji/ LpnujorQ^^/^^ as a 
fimction of r/rp. We show three curves corresponding to the axis-symmetric modes (n = 0) with a fixed value of k 
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(fc = 0.2(27r/ao) and three possible frequencies (w = 174.8, 201.9, 238.1 cm"-'^) involving three different modes (the 
modes LO4, TO3, and LO3 respectively). In Fig. 3 the oscillatory behavior of the interaction potential strength is 
seen, a clearly understandable fact if we realize that it is given as a combination of Bessel functions. At the point 
r = the quantity is not zero, while the points of maxima determine the values of r where the interaction strength 
should be stronger. Moreover, for the modes with increasing values of the frequency the strength of the interaction 
gets weaker. When plotting the absolute value of div (u) we have lost all information concerning its dependence with 
respect to 9 and z. However, we retained the important information about its r dependence. In should also be kept 
in mind that, in order to obtain a better estimation of the electron-phonon interaction strength, it is necessary to 
evaluate the corresponding matrix elements involving the electron wave functions. If the electron wave function is 
peaked at a point where the deformation potential has a low value, we should realize that the interaction strength still 
may be strong depending on the relative weight of both effects. The plot displayed in Fig. 4 is quite similar to that 
of Fig. 3, but now considering the n ~ 1 case (flexural modes). We again see the oscillatory behavior for the same 
reason as in Fig. 3, but now the involved quantity is zero at r = 0, and the calculated strengths are somewhat lower. 
It is also seen that the deformation potential is not lower for higher values of the frequency. In fact, the intermediate 
frequency mode leads to weaker intensities that the lower frequency mode. The same comments as those done for 
Fig. 3 can be made. 

In conclusion, we have calculated the optical phonon modes for a Ge nanowire in the framework of a continuum ap- 
proach. The applied theoretical treatment was adapted to the non polar case, and also was improved by incorporating 
the important linear terms in the wave vector. In this context we expect the theory to provide better results than those 
given by other long wavelength treatments. But, as was already remarked, we do not presently have experimental 
results or microscopic atomistic calculations allowing a desirable comparison with our findings. Nevertheless, and due 
to the current interest in the study of this kind of nanostructures, we think the results of the present work should be 
useful in the way to get a better understanding of Ge and Si nanowires. 



Appendix A. 

The nine functions Gij , appearing in Eq. (|18() , are given in the form 



G23 
G31 

G32 



G: 



33 



2ifc/3| 



Gil — PlQl 
G12 

Gi3 

G21 
G22 



iMv) - vJniv)] , 



7]- 



iPlk [2-a/y + hro/{f3lv)\ J^irj) 



0^ 



[l + 6ro/(/3|j^)] J„(m), 



(Al) 



where a = ^tq, v = fcrp, rj — QlTq, and ^ — Qt'Tq- On the other hand, the functions Si, present in Eqs. H20() . H21() 
and l(22Jl, are defined as 



Si = G32G23 — G22G33 , 

^2 — G21G33 — G31G23 , 
'S'3 = G31G22 — G32G21 . 



(A2) 
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